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The hydrodynamic description of a superfluid is usually based on a two-fluid picture. We compute 
the basic properties of the relativistic two-fluid system from the underlying microscopic physics of 
a relativistic ip 4 complex scalar field theory. We work at nonzero but small temperature and weak 
coupling, and we neglect dissipation. We clarify the relationship between different formulations 
of the two-fluid model, and how they are parameterized in terms of partly redundant current and 
momentum 4-vectors. As an application, we compute the velocities of first and second sound at 
small temperatures and in the presence of a superflow. While our results are of a very general nature, 
we also comment on their interpretation as a step towards the hydrodynamics of the color-flavor 
locked state of quark matter, which, in particular in the presence of kaon condensation, appears to 
be a complicated multi-component fluid. 



I. INTRODUCTION 



A superfluid at finite temperature can be described "macroscopically" in terms of a two-fluid model involving a 
superfluid and a normal fluid. It can also be described "microscopically" in terms of a condensate associated with 
the spontaneous breaking of a U(l) global symmetry in an underlying field theory. In this paper, we show how to 
relate the macroscopic hydrodynamics of a non-dissipative relativistic superfluid to the microscopic field-theoretic 
description. 



A. Macroscopic two-fluid description 

The two-fluid picture was suggested by Tisza [T] shortly after the superfluid nature of 4 He was established experi- 
mentally [2, 3 . At zero temperature, there is only a superfluid component, while at nonzero temperatures below the 
critical temperature, both fluid components coexist. Landau [4] predicted a "second sound" mode where the relative 
densities of the two fluids oscillate, resulting in a temperature (or entropy) wave. This is in contrast to an ordinary 
sound wave, where the total density (or pressure) oscillates. 

The generalization to relativistic superfluids was introduced by Lebedev and Khalatnikov [5j [6] and Carter [7] . The 
approaches in these works — termed "potential" and "convective" variational approaches, respectively — differ in their 
formulation but are equivalent and can be translated into each other [5J 15] . The convective approach uses two four- 
currents, corresponding to total conserved charge and entropy (which is also conserved in the dissipationless case). 
The potential approach uses two conjugate four-momenta, which correspond to the superfluid four-velocity (multiplied 
by the chemical potential) and a thermal four- vector whose temporal component is the temperature. The conventional 
presentation of the two-fluid model, although it looks like two currents (super current and normal current) is really 
based on one current (the entropy current) and one momentum (proportional to the superfluid velocity). We explain 
how the various formalisms are related via Legendre transforms of two Lorentz scalar quantities, a generalized pressure 
and a generalized energy density. We calculate all these quantities in terms of the microscopic theory. 



B. Microscopic field-theoretic description 



The main physical ingredient of a superfluid is a Bose-Einstein condensate, which is responsible for the frictionless 
flow. This was realized by London, who pointed out the connection between the phenomenology of 4 He and the 
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quantum statistics of bosons [TU]. In a more modern, field-theoretic language, the superfluid state spontaneously 
breaks a £7(1) symmetry of the underlying Lagrangian. As a consequence, there is a Goldstone mode whose dispersion 
is linear for small momenta. Using the terminology introduced in the context of 4 He, this mode is called a phonon. 
Roughly speaking, thermal excitations of this mode constitute the normal component in the two-fluid picture, while 
the superfluid component corresponds to the quantum mechanical ground state. We will clarify this distinction below 
(see Sec.[nB|). 

Our microscopic model will be a relativistic ip 4 field theory with a £7(1) symmetry that is spontaneously broken by 
a Bose-Einstein condensate of the fundamental scalars. We derive the two-fluid picture from this model. We express 
the basic currents and momenta of the two-fluid formalism in terms of field-theoretic quantities, and calculate the 
effective action at nonzero temperature in the presence of a homogeneous, static superflow. 



C. Previous work and future applications 

Our approach is related to several existing studies in the literature. For instance, in Ref. jTT] the two-fluid formalism 
is connected with a simple statistical approach to the phonon contribution, Ref. [T2] connects it with a Walecka model 
describing nuclear matter in a neutron star, and Ref. [13] makes the connection to a very general effective field theory. 
In Ref. |14j . a hydrodynamic interpretation of a field-theoretic effective action was discussed for zero temperature. 
For the simplified case of a dissipationless, homogeneous fluid, our study is a generalization of this work to nonzero 
temperatures. In that reference, also an effective Lagrangian for the phonons is formulated which has been employed 
to study transport properties of quark matter [T5H2Ti] . 

Our study is very general since we do not have to specify the system for which our "microscopic" theory, the 
cp 4 model, is an effective description. As a specific application, however, we are ultimately interested in a complete 
hydrodynamic description of the color-flavor locked (CFL) phase of dense quark matter [3T]. The hydrodynamics of 
dense matter, be it dense nuclear or quark matter, is very important for the understanding of the physics of compact 
stars. One prominent example that relates the hydrodynamic properties to an observable - in this case the rotation 
frequency of the star - is the r-mode instability [22] , for recent developments relating the instability to microscopic 
physics see Refs. [2"3H2"o] . Other examples include pulsar glitches, magnetohydrodynamics, and asteroseismology; 
see for instance Refs. 27 Til for applications of multi-fluid hydrodynamics in the astrophysical context. Since the 
densities in the cores of compact stars can be as large as several times nuclear ground state density, it is conceivable 
that the cores are made of deconfined quark matter. The CFL phase is the ground state of three-flavor, asymptotically 
dense quark matter and therefore also of fundamental theoretical interest. Whether it persists down to densities that 
exist inside a compact star is unknown, and many other candidate phases have been discussed 32J. Taking into 
account the effect of a non-negligible strange quark mass, the next phase down in density appears to be the kaon- 
condensed CFL phase (CFL- AT ) [33], whose superfluid properties are expected to be nontrivial. Firstly, the quark 
Cooper pair condensate of CFL spontaneously breaks the £7(1) associated with baryon number conservation, making 
CFL a baryon superfluid. Secondly, the kaon condensate in CFL-A spontaneously breaks the £7(1) associated with 
strangeness conservation. This £7(1), however, is explicitly broken since weak interactions have to be taken into 
account in the astrophysical context. As a consequence, the corresponding Goldstone mode has a mass of the order 
of 50 keV [34] which is small compared to the critical temperature for kaon condensation, estimated to be of the order 
of tens of MeV [35]. Therefore, the question arises how many (super)fluid components CFL- A has, and how the 
two-fluid picture is modified (or whether superfluidity survives in any sense) when there is a small explicit symmetry 
breaking of the £7(1). 



D. Overview 



In Sec. [TT] we start with the zero-temperature case, where the relation between field theory and superfluid hy- 
drodynamics is well established. We elaborate on this relation, which sets the notation and prepares us for the 
nonzero-temperature case. The remaining paper is divided into four parts: the microscopic nonzero-temperature cal- 
culation, Sec. |III| introducing the two-fluid formalism, Sec. |IV[ and connecting both, Sec. |Vj then we use our results 
to compute the sound velocities in the presence of a superflow in Sec. |VI| 

Our convention for the metric tensor is g^ v — diag(l, —1, —1, —1). Our units arc h = c = ks = 1- We work in the 
imaginary-time formalism, i.e., T/VJ2k = -^Sn/ ^ 3 k/(27r) 3 , where n labels the Matsubara frequencies w„ = ik^. 
For bosons oj n = 2nirT, for fermions w„ = (2n + l)irT. 
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II. ZERO TEMPERATURE: SINGLE FLUID 



A. Lagrangian 



Our starting point is the Lagrangian 



C = d fl <pd» ( p*-m 2 \Lp\ 2 -\\<p\\ 



(1) 



with the complex scalar field tp, the mass parameter m > 0, and the coupling constant A > 0. The Lagrangian is 
obviously invariant under global U(l) rotations cp — > e la ip which leads to a conserved charge 1 . Later, in our explicit 
low-temperature results, we shall set m = for simplicity. Note that, since m 2 > 0, the quadratic term of the 
classical potential is positive, and spontaneous symmetry breaking can only occur if we introduce a chemical potential 
associated to the conserved charge. If this chemical potential is larger than m, a Bose-Einstein condensate of the 
particles is formed. In our formalism, we shall introduce such a chemical potential through the time dependence of 
the phase of the order parameter. 

We allow for Bose-Einstein condensation in the usual way by separating the condensate, 



<p{x) 



V2" 



\p{x) + ip\{x) +iip' 2 (x)] . 



(2) 



Here, p(x) is the modulus and ip{x) the phase of the condensate. For convenience, we have introduced the transformed 
fluctuation field <p'(x), which we have written in terms of its real and imaginary parts. Inserting this into the 
Lagrangian yields 



with the tree-level potential 



U 



£=-[/ + £(!) +£(2) +£(3) +£« 



\d*P&P - y - ™ 2 ) + \p* 



and the fluctuation terms, listed by their order in the fluctuation from linear to quartic, 

0,Mpd^2 - ^P) + P(d^d^ - m 2 - A P 2 M , 



£(2) 



(3) 



(4) 



(5a) 



(5b) 



£(3) 



-Xpip[(ipf 



(5c) 



£(4) 



A 



i ,J2\2 



(5d) 



In this paper, we shall only need the tree-level potential and the terms quadratic in the fluctuations because we restrict 
ourselves to the one-loop effective action, see discussion at the beginning of Sec. |III| In the approximation of a static 
and homogeneous superflow that we shall use for our microscopic calculations throughout the paper, the linear terms 
in the fluctuations vanish upon space-time integration in the action because the fluctuations tp' 1: ip' 2 do not contain a 
zero- momentum Fourier component, which is separated by the condensate. 



In the context of CFL-K , this Lagrangian can be read as an effective theory for (neutral) kaons (neglecting the small explicit symmetry 
breaking due to the weak interactions). In this case, m and A can be expressed in terms of the effective kaon chemical potential and the 
effective kaon mass which, in turn, can be related to high-density QCD via the high-density effective theory |33| . 
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B. Zero-temperature hydrodynamics 



Let us first discuss the tree-level potential and ignore all fluctuations. This allows us to present the translation to 
the hydrodynamic equations and quantities in the simplest case, based on Ref. |14j . Later we shall include fluctuations 
to generalize this translation to nonzero temperatures. 

The classical equations of motion are 

Up = pid^d^tp - m 2 - Xp 2 ) , (6a) 

d^p 2 d^j) = 0. (6b) 
The held-theoretic expressions for the conserved charge current and the stress-energy tensor are 

T „„ = 2 i(^£) =2 «£__ 
where we have abbreviated 

a 2 = d^ipd"^ . (8) 

For the stress-energy tensor, we have used the gravitational definition by formally introducing a general metric g 
which, after taking the derivatives, we set to be the metric of flat Minkowski space. This definition guarantees that 
T^ v is symmetric and conserved. Since is also conserved, we have 

d„r = Q, 8^=0, (9) 

which are the hydrodynamic equations. As we expect from Noether's theorem, current conservation is identical to 
the equation of motion (6b) for tp. 

As an ansatz for the solution of the equations of motion we may choose 

V> = p^+Re^^ke i(w *- k ' x) , (10a) 



P = J^^+Re^^"*). (10b) 



In this ansatz, we assume d^rp and p are each composed of a static part plus small oscillations around it. The static 
part of the solution is the superfluid mode, corresponding to an infinite and uniformly flowing supcrfluid. The density 
and flow are specified by the values of the components of p^, which are pure numbers, not functions of x, and are not 
constrained by the equation of motion. The value of p^ is determined by the boundary conditions, which specify the 
topology of the field configuration, namely the number of times the phase winds around as we traverse the space-time 
region in which the superfluid resides. The oscillations around the static solution yield two modes whose dispersion 
relations w(k) are determined by the equations of motion. One of the modes is the massless Goldstone mode, the 
other is a massive "radial" mode. The oscillatory modes can be thermally populated to yield the normal fluid. The 
boundary conditions place no constraint on the amplitudes of the oscillatory modes because they are topologically 
trivial. 

In this section we will assume uniform density and flow of the superfluid, but in general one could obtain a space- 
time-dependent superfluid flow by giving nonzero classical background values to the Sipk and Sp^. The normal fluid 
would then consist of a thermal population on top of that background. 

In Sec. Ill A we will use the poles of the tree-level propagator to determine the dispersion relation of the oscillatory 



modes. This yields the same result as solving the equations of motion ^ with the ansatz (10 1. Oscillations in ip will 



also be relevant to the calculation of the sound velocities in Sec. |VI| In fact, the coefficient of the linear term in the 
massless mode w(k) is nothing but the velocity of first sound. For now, we set 5ipk = Sp^ = 0. In this case, p M = d^tp, 
and we can write 



(11) 
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Inserting this solution back into the Lagrangian gives 



C = -U = 



(a 2 - m 2 f 
4A 



(12) 



It is clear from the usual definition of the pressure in terms of the partition function that, in our current tree-level 
treatment, this Lagrangian is identical to the pressure. 

We now explore the connection between field-theoretic quantities and hydrodynamics. The hydrodynamic approach 
involves writing the current and stress-energy tensor in terms of fluid-mechanical quantities 



t» v = ( £s + p s yv - gTP. , 



(13a) 
(13b) 



where is the superfluid velocity; by definition v^v^ = 1. To identify n s , e s , and P s we go to the superfluid rest 
frame = (1, 0, 0, 0), where there is no flow of charge, energy, or momentum: j l = T° l = 0. Then we have j = n s , 
T 00 = e s , T lJ = 5 1: >P S , so n s , e s , and P s are the charge density, energy density, and pressure in the superfluid rest 
frame; the pressure is isotropic. (We will discuss the superfluid density p s below). They can be expressed covariantly 
via appropriate contractions, 



6 S = VaVvT^ , 



(14) 



To relate these hydrodynamical quantities to the microscopic physics we can use Eqs. ([7]). For n s and w M we obtain 
immediately from ( 7a ) , 



A 



7(!.v s ) , 



(15) 



where v s is the superfluid 3- velocity with Lorentz factor 7 = l/yl — vf. The velocity is proportional to the gradient 
of a scalar field, v s = —^ip/doip, which is what one expects since a superfluid without vortices is irrotational. [The 
minus sign appears because the 3- velocity v s corresponds to the spatial components of the contravariant 4- vector v^, 
while the operator V corresponds to the spatial components of the covariant 4-vector d^, i.e., = (dt, —V).] 
For e s and P s we first use (14) and (7b) to obtain partly microscopic expressions 



Ps 



C + (v^ib - a)n s 



(16) 



From contracting the relation for in (15) with d^ip we know that v^d^ip — °~\ as a consequence, we obtain the 
expected relation P s — C. We can now identify the physical meaning of a: using the zero-temperature thermodynamic 
relation in the superfluid rest frame e s + P s = /i s n s , we find 



fi, 



= <T = V^dnlp , 



(17) 



so a is identified with fi s , the chemical potential in the superfluid rest frame. Going back to the solution for the 
modulus of the condensate (11) we see that, as expected, Bose condensation only occurs for fi 2 > m 2 . Using (12) and 
( 16 ) we finally obtain fully microscopic expressions for the energy density and the pressure, 

(a 2 - m 2 ) 2 



P, 



4A 



(3cr 2 + to 2 ) (a 2 -m 2 ) 
4A 



(18) 



3P S 



2 2 

TO p 



Note that to is the only mass scale in our Lagrangian. The trace of the stress-energy tensor = 
vanishes for to = 0. 

Finally, the superfluid density p s is defined via the expansion in small three- velocities of the momentum and energy 
densities 



nOO 



T 0t = Ps v S! +O(|v s | J ), T u 
From these expansions we obtain its microscopic form 

p s = e s + P s 



PsVi 



o(Ws\ 4 ) 



(19) 



(20) 
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C. Single-fluid formalism 



Let us now, following Refs. [SHS], introduce the formalism that we shall use below to describe the two fluid com- 
ponents at nonzero temperature. Although this might seem unnecessary for the zero-temperature case where there is 
only a single fluid, it is a useful preparation for the nonzero temperature case. 

In the hydrodynamic formalism, the basic variables that describe one fluid are j M , which is the conserved current 
density, and d^ip, which we will identify as the conjugate momentum. Motivated by Eqs. Q, we write the stress-energy 
tensor as a function of these variables 



JVM/ = _ g »U^ + jHQK^ ^ (21) 

where we call W the generalized pressure, and we define the generalized energy density 

A = + 3* = + fd^ . (22) 

By comparison with the expressions from the previous subsection it is easy to see that in the single-fluid case ^ and 
A are simply pressure and energy density in the fluid rest frame, = P s = £, A = e s . It is then natural to treat A 
and ^ as Legendre transforms of each other: 'J depends only on d^ip, and we assume that A depends only on (the 
fact that OA/did^) = follows from (22]) and ((7a|). So 

dA = df , = fd(d^) . (23) 

Since A and VP are Lorentz scalars, they depend on the Lorentz scalars j 2 = j^j^ = n 2 s and a" 2 — d^ipd^ip, respectively. 
This allows us to write 

fi^ = |£ = *j/.. ^ 2 f^ (24) 

and 

r = am =Bffl *> (25) 

These relations say that the current is proportional to the conjugate momentum with a coefficient given by the 
underlying microscopic physics. In this one-fluid system, we simply read off 

B = B- 1 = g2 ~ m2 . (26) 
A 

This formalism will allows us to introduce the so-called entrainment coefficients as the coefficients in a two-fluid 
system which relate the current of the first fluid with the conjugate momentum of the second fluid and vice versa. 
But before we do so, we go back to our microscopic calculation to "prepare" it for the translation into the two-fluid 
picture. 



III. NONZERO TEMPERATURE: MICROSCOPIC CALCULATION 



A. Effective action and anisotropic dispersions 



In our context, the nonzero-temperature calculation is complicated for two reasons: (i) we are interested in the 
situation of a nonzero (homogeneous) superflow; this makes the system and in particular the dispersions of the 
Goldstone mode anisotropic and (m) in order to compute the properties of the system for all temperatures up to the 
critical temperature, a self-consistent formalism such as the Cornwall- Jackiw-Tomboulis (CJT) formalism [5BH55] is 
needed. (For a detailed discussion of this formalism in the same context, but without superflow, see Ref. |35j.) 

Here we are mostly interested in analytical results to discuss the translation to hydrodynamics thoroughly. There- 
fore, we shall restrict ourselves to small temperatures and ignore the complications from point (ii). This means that 
we shall not determine the condensate self-consistently, and rather work with its zero-temperature value. Starting 
from the full self-consistent formalism, one can show that the temperature-dependence of the condensate includes 
an additional power of the coupling constant compared to the terms we are keeping. We are thus working in the 
low-temperature, weak-coupling approximation. Note, however, that the melting of the condensate, as well as other 
contributions from the full action may induce terms proportional to \/j, 2 T 2 in the pressure, while we shall only keep 
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terms proportional to T 4 and higher order in T. This approximation is only consistent if T 2 ^> A// 2 , i.e., strictly 
speaking, our approximation leaves a "gap" between T = and the small temperatures we are discussing, although 
this gap can be made arbitrarily small by choosing the coupling constant A sufficiently small. We defer a careful 
discussion of the self-consistent CJT calculation with superflow, including a numerical evaluation up to the critical 
temperature to a forthcoming publication |39j . 

Within this approximation, we can work with the simple one-loop effective action 



r 



V 
T 



U(p,a) 



y2 



(27) 



with the tree-level potential U from Eq. Q and the inverse tree-level propagator in momentum space S~ 1 (k) which 
can be read off from the terms quadratic in the fluctuations in Eq. (5b). The trace is taken over the 2x2 space of 
the two real degrees of freedom of the complex scalar field, and the sum is taken over four- momenta fc M = (fc ,k), 
ko = —iu) n with the bosonic Matsubara frequencies u n . The sum over three-momenta, here written for a finite volume 
V and thus over discrete momenta, becomes an integral in the thermodynamic limit. At the stationary point, i.e., 
with Eq. ( 11 ), we have 



S~\k) 



-k 2 + 2{<r 2 - to 2 ) 2ik-dip 
-2ik ■ dip -k 2 



(28) 



Here, k = k^k^, k ■ dip = k^d^ip and a given in Eq. ([8|). In all our microscopic calculations we evaluate thermal 
fluctuations in the background of a uniform superflow, so we take d^ip (and hence a) to be space-time independent. 
Since we neglect the melting of the condensate, we have inserted the zero-temperature solution into the propagator. 

In order to compute the effective action, one may employ partial integration with respect to the |k| integral. The 
effective action density then becomes 



T 



2^2 



(cr 2 — TO 2 ) 
4A 



2T^k 2 (fc 2 

fc 



3V 



a 2 + to 2 ) + 2k- dip k-Vip 
det 5-!(fc) 



(29) 



where we have used the explicit form of the tree-level propagator and have performed the 2x2 trace. We have also 
inserted the solution (11) into the potential U. The evaluation of the Matsubara sum and the simplification of the 
result is explained in detail in appendix [X] Neglecting the thermal contribution of the massive mode, we can write 
the result as 



(a 2 — m 2 



d 3 k 



^(ei.k.k) 



V 4A J (2tt) 3 (ei, k + ei,_ k )(ei. k + e 2 ,-k)(ei,k-e2,k; 

where -F(fco, k) denotes the numerator in the momentum sum of Eq. ( |29[ ), 

F(k ,k) = - 2 [k 2 (k 2 - a 2 +m 



coth 



£i,k 
2T 



-2k ■ dipk-Vip] , 



and where £i/2.k & re complicated excitation energies whose small- momentum approximations are 



ei.k 



3fT 2 



C(k)|k| + o(|k| 3 ) 



e 2 ,k = V2^Sa 2 - to 2 + 2(W) 2 + 0(|k|) . 



Here we have abbreviated 



C(k) = 




(Vip) 2 - (Vtp ■ k) 2 

3<7 2 — TO 2 



2d ipVip ■ k 



Va 2 ^ 



TO 2 V / 3(T 2 



2(vyp 2 

3<7 2 — TO 



The physically relevant low-energy excitation is e^k- This is the Goldstone mode which, as Eq. (32a 
massless and linear in the momentum for small momenta, as it should be. The coefficient in front of the 
the speed of (first) sound. For the case without superflow we have C(k) = 1, and thus, in the limit to = 



(30) 



(31) 



(32a) 



(32b) 



(33) 

) confirms, is 
linear term is 
0, we recover 



the well-known The superflow introduces an angular dependence C(k) into the Goldstone dispersion and thus 



also into the sound velocity. This angular-dependent function also shows that, in contrast to the zero-temperature 
case, we cannot write the result in a covariant way since the temporal and spatial components of d^ip appear separately. 
Another observation is the complicated factor in the integrand of Eq. ( p0| ) . This factor indicates a mixing between the 
original modes of the complex field (between particles and antiparticles, essentially) due to condensation. Such a factor 
appears also, although considerably simpler, for the case without superflow. It is analogous to a Bogoliubov coefficient 
in the case of Cooper pairing which, in that case, accounts for the mixing between fermions and fermion-holes. 
We shall further evaluate the effective action for small temperatures in Sec. |III C| 



B. Stress-energy tensor and current 



The stress-energy tensor and the current are given by 

2 &{y/=gC) 



J 



DC 



(34) 



where the angular brackets denote the expectation value in the finite-temperature ensemble of the microscopic theory, 



(A) = iy VLp'{Dif' 2 A exp (J d 4 x^gC^j , 



(35) 



with the partition function Z defined so that (1) = 1. After performing the functional integral we can write the results 
as 



dg^ J V ^ 



S(k) 



dS^jk) _ g^_ 
dg^u 2 



(36a) 



A 



1 T 
2V 



5> 



S(k) 



aS" 1 (A:) 



(36b) 



In the case of the stress-energy tensor, one can check explicitly that the sum over the Matsubara frequencies leads to 
an infinite result. A renormalization is thus required. As a renormalization condition, we require for the case without 
superflow the obvious interpretation of the diagonal components T^ v in terms of the energy density e and the pressure 
P, 



T 



no 



T l ° = 8 lJ P 



(37) 



Then, switching on a nonzero superflow, does not yield any additional divergences. These conditions can be imple- 
mented on a very general level, without explicit evaluation of the stress-energy tensor. We do so in detail in appendix 
|B] The calculation in this appendix also leads to a very useful formulation of the effective action, the stress-energy 
tensor, and the current, which will later facilitate the interpretation in terms of the hydrodynamic two-fluid picture. 
We define (for motivation see Eq. (72)) 



* fe = * fe (cr 2 , k ■ dip, k 2 ) ee — Trln 



s-\k) 



J-2 



and 



d(k-d^p) detS- 1 



(38) 



(39a) 



Bk = 2 



C k 



d^ k _ 2k 2 
~~ detS- 1 



' ' dk 2 



2(k 2 



det S~ 



(39b) 



(39c) 



The reason to introduce these quantities and the choice of notation will become clear after we have discussed the 
two-fluid formalism. In the zero-temperature discussion of the single-fluid dynamics we have already encountered the 
generalized pressure "J, of which the notation is reminiscent. The notations Ak, Bk, Ck are chosen in analogy to 
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A, B, C, to be introduced in Sec. IV A (we already know the coefficient B from the single-fluid treatment in Sec. II C I 



With the help of these quantities we can write the action as (see appendix [BJ) 



V 



(40) 



where we have abbreviated it' 1 = (1, 0, 0, 0), while the renormalized stress-energy tensor and the current become 



V" = ( 2 - //' J 7 

og^v 



(41a) 



3 



0*V 



g - m 1 T 



V 



X 



(B k d»1> + A k ¥ 



(41b) 



It is instructive to compute the trace of the stress-energy tensor. With the help of Eq. (|B3j) we immediately find 

A- 2 



T"„ = m 2 



k 



(42) 



i.e., the trace vanishes when we set m = 0, as expected. 



C. Evaluation for small temperatures 



We want to compute the effective action, the components of the stress-energy tensor, and the current for small 
temperatures explicitly. To this end, we have to approximate momentum sums of the form 



T F (k) 
V ^ detS-^fc) ' 



(43) 



where, for the case of the effective action, F(k) is given by Eq. (31 ), and for the stress-energy tensor and the current 
we need to replace F(k) by 

F^(k) = 2[-{k 2 -a 2 +m 2 )k^k v + k 2 d^^d v ^ + 2{k-d^){k^d v ij + k u d^i!) + u^u v detS- 1 {k)\ , (44a) 



F»(k) = 2{k 2 d^ip + 2k ■ dip . 



(44b) 



which can be read off from Eqs. (41) 



For all three cases, we can write the result of the Matsubara sum in the form (30 1. Then, we write coth[ei i k/(2T 1 )] 



1 + 2f{e\^) with the Bose distribution function f[x) = \j(e x l T — 1). The integral over the first term is, in our 
approximation, temperature independent. It is divergent and has to be renormalized by subtracting the vacuum 
contribution. After doing so, a finite term remains which however is suppressed by one power of A compared to 
the zero-temperature term oc A" 1 we have already computed, for instance in Eq. (30). We shall thus neglect this 



contribution and only keep the thermal contribution, i.e., the integral over the second term that contains the Bose 
function. This contribution is finite and unaffected by the renormalization. 2 

We explain the small-temperature expansion in detail in appendix [Cl Before we come to the main results, let us 



discuss the simpler example without superflow, VV> = 0. In this case, e^k = ejj", €2. k = c k with 



\Jk 2 + 3^i 2 - m 2 =F \/^\i 2 k 2 + (3/i 2 - m 2 f , 



(45) 



2 The sublcading zero-temperature contribution we are neglecting here gives rise to the difference between the superfluid density and 
the condensate density. At zero temperature, the charge density is always identical to the superfluid density, n(T = 0) = n s , while 
the condensate density, defined as fip 2 with the modulus of the condensate p, might be smaller. In our leading-order approximation, 
condensate and superfluid densities are identical. 
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,,4 

£t< 1 -'» 


ir 2 T 4 1 If 2 


47T _Z 1 — V s 


10^3 (1 - 3v?)3 


105^ M 2 (l-3v?) 6 




1 2 

1 - v s 


(1-v2)(1-3v2) 


-(1 - v 2 )(l - 3v 2 )(5 + 30v 2 + 9v 4 ) 


j,00 


3 + v 2 


3 - 20v 2 + 9v 4 


-(15 - 160v 2 - 774v 4 + 432vf + 135vf ) 


rpOi 


4v sl 


-8v si 


2v sl (95 + 243v 2 - 135v 4 - 27v!j) 


T ± 


1 2 

1 - v s 


(l-v 2 )(l-3v 2 ) 


-(1 - v 2 )(l - 3v 2 )(5 + 30v 2 + 9v 4 ) 


T|| 


l + 3v 2 


1 - 12v 2 + 3v 4 


-(5 - 180v 2 - 582v 4 + 324vf + 81vf) 


«° 


4 


-8v 2 


2(5 + 105v 2 + 147v 4 - 81vf ) 




4v s 


-8v s 


2v s (95 + 243v 2 - 135v 4 - 27vf) 



TABLE I: Microscopic results for the effective action V, the stress-energy tensor T 1 "*, and the conserved charge current 
up to order T 6 . Each row is the result for the quantity given in the left column; this quantity is a sum of the (i 4 , T 4 , and 
T 6 /u 2 terms given in the top row, each multiplied by the specific entry of the table. For example, the first row is equivalent 
to Eq. ( |47| . These results are obtained without making any assumptions about the magnitude of the superfluid velocity v s . 
We see that there is a divergence in all nonzero temperature results for v 2 — > 1/3, indicating an instability of the superflow, in 
accordance with Landau's critical velocity for superfluidity. 



where we have set d^tp = \i. Here, t\ is the Goldstone mode, while e k is the massive mode. This yields the 
small-temperature result for the pressure to order T 6 (see appendix [c]) , 

V { V ' ~ 4A ( M 2 -m 2 ) 3 / 2 90 (^ 2 -m 2 ) 7 /2 63/i 2 ' 1 ' 

The T 4 term has two interesting properties related to the resulting charge density, which is obtained by taking the 
derivative with respect to fi. Firstly, if we set m — 0, the /x-dependence drops out, such that there is no T 4 //i 
contribution to the charge density in this case. Secondly, in the presence of a finite m, one finds that the T 4 //i 
contribution to the charge density is negative, i.e., for small temperatures the density decreases with temperature. 
This somewhat surprising behavior c an b e traced back to the dependence of the sound velocity on the chemical 



potential: the T term in the pressure (46) can be written as it T /(90wj), where u\ = (f/r — m ) ' /(3p — m ) ' is 



the speed of first sound which increases monotonically from zero for /i = m to the asymptotic value l/v3 for /j, — > oo. 
This reflects the fact that as the speed of sound increases, the range of momenta corresponding to phonon energies of 
order T gets smaller, so the number of phonons in the normal fluid goes down, so the T 4 term in the pressure decreases 
monotonically with /i, and thus the charge density is negative for all fi > m and grows more negative with increasing 
temperature. (The second derivative d 2 P/dfj? is positive, i.e., the system is thermodynamically stable.) The T 6 term 
has neither of these properties, it contributes to the charge density even for m = and gives rise to a positive T 6 //it 3 
term in the density. The case with superflow is of course more complicated. In particular, the momentum integration 
now involves a nontrivial angular integral over the angle between the momentum and VV>- Nevertheless, it turns out 
that this angular integral can be performed analytically for all cases we consider. For brevity, we set m — in the 
following. Moreover, we denote jx = d^ip. At this point, this may be viewed as a mere notation; we shall see below 



that /i is the chemical potential in the normal-fluid rest frame, which is defined below Eqs. (49). Then, we obtain for 
the effective action density to order T 6 

T T~»\l v 2 ! 2 I (1 ' V ' )2 ^ (l-v 2 ) 2 f5| oq^.q^ (m 

v r -AX {1 -^> + Trvl (i - 3v 2 ) 2 (i - sv 2 )* (5 + 30v * + 9vJ • (47) 



Obviously, for v s = we recover the m = limit of Eq. ( 46 ) . Note that the result is valid to all orders in the superfluid 
velocity. We have not applied any expansion in |v s |. Analogously, we compute the components of the stress-energy 
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tensor and the current. For the spatial components of the stress-energy tensor we define 



T 



(vv-) 2 



T, 



11 - (V^) 2 ij 



(48) 



which the spatially transverse and longitudinal components with respect to v s . We collect all results in Table [T] 

All expressions in the table are written in terms of frame-dependent quantities. At zero temperature it is very 
natural to write them in terms of the Lorentz scalar a = /Ltyl — v 2 , for instance T± = c 4 /(4A), as we have done 
in Sec. [TTJ Such a formulation is less obvious for the nonzero temperature terms. As we shall discuss in Sec. [Vj the 
microscopic calculation has a preferred rest frame which we shall identify as the normal-fluid rest frame, and all 
thermodynamic variables in the above expressions, i.e., /i, T, and v s , are measured in this particular frame. As a 
consequence, the velocity-dependence shown here is a mixture of trivial Lorentz boosts and complicated effects of the 
superflow on the collective excitations. 



IV. RELATIVISTIC TWO-FLUID FORMALISM 



So far we have discussed the zero-temperature limit of the microscopic theory and its formulation in hydrodynamic 
terms, as well as the generalization to nonzero temperatures of our field theory. Before translating the latter into 
the language of hydrodynamics we need to introduce the two-fluid formalism. As indicated in the introduction, a 
superfluid becomes a two-fluid system for temperatures larger than zero and smaller than the critical temperature. At 
zero temperature there is only the superfluid component, while at the critical temperature (and above) there is only 
the so-called normal component. The discussion in this entire section is general in the sense that we do not assume 
a uniform superflow; we shall get back to this assumption only in the next section where the microscopic results are 
needed. The most straightforward way to add a second component to the single-fluid hydrodynamics seems to be 



n n u^ + n, 



(49a) 



= (e. 



P n )uV - gTP n + (e s + P.)^^- _ g ^ Ps 



(49b) 



Here we have simply taken the zero-temperature expression ( 13 1 with the superfluid velocity ( 15 1 and added analogous 
terms for the normal component. In particular, we have introduced the normal fluid velocity u M , with u^u^ = 1. As 
a consequence, we can now also define the normal-fluid rest frame by = (1, 0, 0, 0). Obviously, if we go to either 
the superfluid or the normal-fluid rest frame, there will be a nonvanishing three-current. In other words, in contrast 
to the zero-temperature (i.e., single fluid) case, there is now in general no frame where the pressure is isotropic. The 
decomposition in terms of superfluid and normal components in the form (49) can be found for instance in Refs. [4"0T - 



|4"4"]. 3 In order to make a first connection to field theory, let us go to the normal-fluid rest frame (as we shall argue 
later, this is the frame in which our microscopic calculation is performed). In this frame, take for instance the spatial 



components of the current (49a), and contract the equation with Vtp. This yields 



(vv>) 2 



(50) 



The right-hand side is now defined in terms of the field-theoretic quantities. We have thus found a microscopic 
expression for the superfluid number density n s . Note that n s is the superfluid number density measured in the 



superfluid rest frame; from the temporal component of Eq. (49a) we see that the superfluid number density measured 
in the normal-fluid rest frame is fin s /a, with \i — d ?/; being the chemical potential in the normal-fluid rest frame. As 
we know from Eq. (15), the factor fi/a = \j \f\ — v 2 is a usual Lorentz factor. 



There is a slight difference in Refs. 1411 1421 . where the stress-energy tensor is written as (adapted to our notation) 



= (fin + Ts)u^u» - g^P + an s 

where n = n„ + ^n B is the charge density, measured in the normal-fluid rest frame, and P = P s + P n . To obtain this stress energy 
tensor from Eq. |49bl, one defines w 11 = d^ip — fiu 11 and uses e s + P s = <rn s , e„ + P n = fin n + Ts. 
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In the two-fluid formulation (49) the charge current is decomposed into two non-conserved currents, so that each 
fluid component is characterized by a non-conserved current. This contrasts with the formulation of Refs. [SHZ], where 
the two components of the fluid are given by the conserved current j M and the entropy current s" (the latter being 
conserved in the dissipationless case). We shall explain this formalism now. This is a generalization of the one-fluid 
discussion of Sec. Ill CI 



A. Two currents, entrainment, and hydrodynamic equations 



Generalizing the zero-temperature, single- fluid stress-energy tensor (21 1, we write 



(51) 



where, as before, j" is the total conserved current (now including thermal contributions), and d^tp the corresponding 
conjugate momentum which, as we know, is related to the superfluid velocity. Analogously, we have introduced the 
entropy current s M and the corresponding conjugate momentum M . The entropy current is related to the four- velocity 
of the normal fluid via 



(52) 



where s = (s^s^) 1 / 2 . 



In the normal-fluid rest frame, s M = (s , 0, 0, 0), so s = s° is the entropy density in the normal- 
fluid rest frame. (Of course, in the superfluid rest frame, s° ^ s). Remember that the generalized pressure \& was, 
in the single-fluid case, identical to the pressure in the superfluid rest frame. Now we may relate ^Sf to the pressures 
of the normal and superfluid components as follows. From Eq. (51) we can derive the following expression for the 



generalized pressure (without using any properties of the various four-vectors) 



dtp(s ll d u tp + s u d^ip) ~ s 2 d ll tpd v tp 
(s ■ dtp) 2 — s 2 a 2 



a 2 s. 



(53) 



Inserting the stress-energy tensor (49b) on the right-hand side of this equation, we find = P s + P n . Consequently, 



the generalized pressure is the sum of the pressures of the superfluid and normal components, each measured in their 
respective rest frames. 

According to its definition A = 



3^ (22), the generalized energy density is now 
A = + j ■ dtp + s ■ 6 , 



(54) 



which can be read as a generalized, covariant thermodynamic relation between pressure and energy d ensi ty. It follows 
that in the normal fluid rest frame 0° is the temperature, although because of entrainment (see ( [56] ) ) the spatial 
components 0* may not vanish. Analogously to the pressure, one can show that A = e s + e n . As already discussed 
in the zero-temperature case, Eq. (54) is a Legendre transform between A and VP that changes the dependence from 
the (now two) currents to the conjugate momenta. The energy density is a function of the currents j", s M , while the 
pressure is a function of the momenta, 



dA = d^dj" + 



(55) 



Since A and ^ are Lorentz scalars, they must depend on the Lorentz scalars built from the currents, s 2 , j , s ■ j, and 
the momenta, a 2 , 2 , O • dtp, respectively. Therefore, we can write 



d^tP 



OA 

dA 
ds,. 



Bf + As" , 



Af + Cs" , 



(56a) 



(56b) 



where 



A 



dA 
d(j ■ s) 



B : 



( <9A 



C 



, dA 
''ds 2 



(57) 



in the notation introduced in Ref. [IT] (where C, B, A are called caloric, bulk, and anomaly coefficients). In general, 
the momenta are not four-parallel to their corresponding currents, but also receive a contribution from the other 
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current. This effect manifests itself in the coefficient A which therefore is also called entrainment coefficient, see for 
instance Ref. [45] . 

For completeness, and because it will turn out to be useful for the discussion in the subsequent section, we also 
introduce a notation for the inverse transformation, 



f = 



= Bd"il> + AW 



d& L 



= Ad^ip + C@L 



(58a) 



(58b) 



with 



A 



B : 



<9(9 • dip) ' ~ da 2 ' 

Obviously, the coefficients are related by a simple matrix inversion, 

B C 



c 



BC-A 2 



B = 



BC-A 2 



A 



A 



BC-A 2 



(59) 



(60) 



In summary, we have written the stress-energy tensor (51 1 in terms of the Lorentz-scalar 'J and 4 four- vectors, two 
of which are independent. Usually, one chooses either the two four-currents or the two four- momenta as independent 
variables. The formulation in terms of normal and superfluid components ( |49[ ) is a "mixed" form, which uses one 
of the currents and one of the momenta as its basic variables, namely s M (which corresponds to u^) and d^ip. To 
translate between these two formulations (see also appendix A in Ref. [43]), we rewrite the current and stress-energy 
tensor with the help of Eq. (56 1, 



(61a) 



BC-A 2 
B 



(61b) 



[As an aside, from this form we see that the stress-energy tensor is symmetric, which is not obvious from the form 
(51 1.] Comparing with Eqs. (49), we can identify 



o 
B 



As 
~~B 



P, 



~B 



P-n 



BC-A 2 
B 



(62) 



With the stress-energy tensor (51) we can also rewrite the hydrodynamical conservation equations We find 







0"VV + j^d^t/j - d v d' J "ip) 



(63) 



where we have used d v ^ = j^d^d^ip + s^d"®^ [because of Eq. (55)]. The first term on the right-hand side vanishes 
due to the current conservation. In our particular case of a superfluid, the second term also vanishes (up to this point, 
the specific form of the momentum d^tp has not been used, i.e., the result holds for arbitrary two-fluid systems). 
As a consequence, at zero temperature where there is no entropy current, the energy-momentum conservation is 
automatically fulfilled due to the conserved, curl- free current. Contracting the remaining two terms with s v shows 
that, due to the antisymmetry of the so-called vorticity d^Q v — d u Q^, the entropy current is also conserved (provided 
that s v Q v ^ 0). Consequently, the conservation equations are equivalent to 



9^ = 0, 



0. 



a^Q" -d v W) = 0. 



(64) 



4 The equation for the vorticity is sometimes (see for instance Ref. jllj ) written with the help of the inverse temperature four- vector 

■ = -— . 
s6 

One can easily check that Sf 1 {d l 'Q^ — 9^0") = s ■ 0(/3 (J c' M © 1 ' + &n,d v /J**), and thus the the vorticity equation becomes 

Pud^e" + e^d"^ = o. 

For a physical interpretation of the vorticity equation, see for instance Sec. 6 in Ref. 9 . 
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normal-fluid rest frame 


superfluid rest frame 


charge density j° 


aV 

rin + n a 

a 


s° 
s 


spatial current j 


n s v s 

(J 


S ° 

-n»v„ 


energy density T 00 


A — J ■ Vlp 


7V -|- S • d 


momentum density 


3 o W 




long, pressure T\ | 


^ — J • v yj 


-[- s * 


transv. pressure T± 






T 00 +T X -T|| 


A 


A 



TABLE II: Components of the current and the stress-energy tensor in the normal and superfluid rest frames. In each frame, 
d ip is the chemical potential, s° the entropy, and 0° the temperature, while n n and n s are the normal and superfluid 
number densities, measured in their respective rest frames and v n = s/s° (v s = ~\7tp/d°4>) the three-velocities of the normal 
(superfluid) component, measured in the superfluid (normal) rest frame. Longitudinal and transverse pressures are defined 
with respect to the three-direction of the velocity of the other fluid component. 



It is instructive to write down the current and the stress-energy tensor in two preferred rest frames, the normal-fluid 
rest frame, defined by = (s°, 0,0,0), and the superfluid rest frame, defined by d^ip = (c^Vj 0; 0, 0). The various 
components are listed in Table [TT| and can be interpreted as follows. 

The components of the current are written in terms of n n and n s . (Alternatively, due to the translation given in 



Eq. (62), we could have written them in terms of the coefficients A and £>.) Since n n and n s are measured in their 
respective rest frames, the charge density of the other fluid component contains an explicit Lorentz factor, i.e., n s is 
multiplied by d°ip/<7 in the normal-fluid rest frame, and n n is multiplied by s°/s in the superfluid rest frame. The 
spatial components of the currents are given by the respective number densities times the three- velocities of the other 
fluid. 

The components of the stress-energy tensor are written in terms of the Lorentz scalars and A. The last two lines of 
the table illustrate the meaning of these quantities. The transverse pressure, i.e., the pressure measured in the spatially 
orthogonal direction with respect to the fluid velocity of the other current, is identical to the generalized pressure 
"J. The energy density T 00 contains the kinetic energy from the other fluid component. This is exactly the term 
that distinguishes the transverse from the longitudinal pressure. Therefore, the combination of the frame-dependent 
quantities T 00 +T± — Tj| is identical to the generalized energy density A. 

We may use this table for a comparison with our microscopic low-temperature results from Table [Tj One can for 
instance check from the explicit results, that the difference between longitudinal and transverse pressure is indeed 
— j • V?/^ = /ij • v s , and that the momentum density T l ° is indeed d°ip — [i times the current j\ More importantly, we 
see that T± = ^T, which, since T± = vp, suggests that the generalized pressure should be identified with the effective 
action. We shall confirm this on a more general level in the next section. 

Finally, let us use the relations summarized in the table to define superfluid and normal densities p s and p n . These 
energy densities are defined in the low- velocity limit, in generalization of the mass densities in the nonrelativistic 
framework (and in contrast to n s , n n which are number densities) . To define the superfluid density p s we first have 
to go into the normal-fluid rest frame. With T {h = j l d°ip, j = — n s V^/a, and v s = —Vip/d°ip, we obtain 

T o* = (dV) 2 ^v 5 , = Ps v sl + 0(|v s | 3 ) , (65) 



with the Lorentz scalar 



p s = an s . (66) 
(Alternatively, we can write p s = a 2 /B.) This expression for the superfluid density is obviously in agreement with the 



zero-temperature result ( 20 1 . The superfluid density appears also in the energy density as part of the kinetic energy, 
T 00 = A + p s v 2 + C(|v s p). The normal fluid density is defined analogously: in the superfluid rest frame, we have 
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T l ° = j l d°ijj + s l Q°. Inserting j l = n n s l /s and using v„ = s/s°, we obtain 



V"= ( ( ;V^H ., ,l B l, )v,„= / -„v,„+CO(|v„| , ) . «»7) 



with 



p n = an n + sQ . (68) 



(Alternatively, p n = — asA/B + sO.) 



V. TWO FLUIDS IN TERMS OF MICROSCOPIC VARIABLES 

How is the two-fluid formalism explained in the previous subsection related to the microscopic calculation of Sec. |III| ? 
In other words, what is the field-theoretic definition of the basic quantities of the two-fluid formalism? We have seen 



in Sec. |II C that for zero temperature the field-theoretic quantities could be translated straightforwardly into (single- 
fluid) hydrodynamics. For instance, the generalized pressure was simply given by the Lagrangian. In particular, the 
microscopic formulation was covariant, even in the presence of a chemical potential. The microscopic results from 
Sec. |III| however, show that this covariance is lost when we go to nonzero temperatures, at least in the explicit results, 
see for instance the dispersion relations from Eq. (32 1. Nevertheless, we shall now show that, at least in our static, 
homogeneous, dissipationless scenario, we may map the field-theoretic formulation onto the hydrodynamic one. Our 
covariant formulation from Sec. |III B| before performing the Matsubara sum will turn out to be very useful for this 
purpose. 



A. Generalized pressure corresponds to effective action 



The first important step is to realize that the microscopic calculation is, in the terminology of the previous subsec- 
tion, performed in the normal-fluid rest frame u M = (1, 0, 0, 0), which is the rest frame of the heat bath. For discussions 
how a general frame can be introduced in relativistic thermodynamics and thermal field theory in a covariant way 
see for instance Refs. [13 U?]. We are not interested in reintroducing a general in our calculation, which would be 
complicated conceptually because we are using the imaginary time formalism and also complicated technically since 
there is already one four-vector in our calculation, namely d^tp. 

Our first goal is to identify the generalized pressure "J '. With the stress-energy tensor (51 1 and using — / s we 
can write 



1 



T 



(69) 



Let us now compute the right-hand side with the microscopic expression for T^" . With T MI/ and j 1 * from Eqs. (41) 
we compute 



\{gT ~ u»u v ){j^ - V) =-U- \(<r - u^u")^ Y, {Ckk^K + A k k^) . 



(70) 



Here we have used — (1, 0, 0, 0), while Eq. ( 69 ) is a general relation for arbitrary four- velocities u M . It is important 
to keep in mind that we cannot simply promote to an arbitrary four-velocity in the microscopic calculation. It 
would occur additionally in different places in the calculation, which we have not identified in our present treatment. 
But, of course, since ^ is a Lorentz scalar, the normal- fluid rest frame is as good as any other frame to compute 
By comparing with Eq. (40) we see that the right-hand side of Eq. (70) is exactly (T/V times) the effective potential. 
Therefore, we have derived our first important result, 



T 

* = — r. 

V 



(71) 



This relation is somewhat expected since we already know that, at zero temperature (and at tree- level), \& corresponds 
to the Lagrangian, which in this case gives the microscopic pressure. At nonzero temperature (without superflow), the 
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A 



B 



d°ip v1j°d°i> + j • Vi/> + vls°e° 



j°sV(v?i aV + j • vvo - j • v^°e° 



(9V) 2 



A 



B 



fgV(y|/gV + j ■ - j ■ v^ s °e° 



TABLE III: Coefficients that relate the currents j M , s M with the conjugate momenta d^ip, O' 1 , given in terms of "microscopic" 
quantities: the Noether current j M , the space-time derivative of the phase of the condensate d^ip, the temperature O , and the 
entropy density s°, all measured in the normal-fluid rest frame. The low-temperature approximations for A, B, C are given in 
Eqs. (761. 



effective action gives the (isotropic) pressure. Therefore, the relation (71 1 is a natural generalization to the anisotropic 
case with a nonzero superflow. This motivates our choice of notation in Eq. ( 38 ) because now we have 



<i- = -U - 



T 
V 



(72) 



Next, let us discuss the generalized thermodynamic relation A = — + j ■ dtp + s ■ 0. In the normal-fluid rest frame, 
s m — ( s o^ o) and thus s ■ = ,s°0 , which is the product of entropy and temperature, measured in this particular 
frame. To confirm this microscopically, we use the thermodynamical definition of the entropy density, 

y E (* fe + 2 + Ckk ° + A k k od ip) , 



df 



(73) 



with A k and C k defined in Eqs. (39). Note that the first two terms and 2) come from the explicit T-dependence in 
the prefactor T/V and in the 1/T 2 within the logarithm. On the other hand, we can compute s ■ 6 via the generalized 
thermodynamic relation. We find 



A + * - j ■ dtp 

T 
V 



E 



2m 2 k 2 
detS- 



*fc + C k k 2 + A k k d ip - (C k k A + B k a z + 2A k k ■ dtp) + 



(74) 



where we have used A = + 3^ (22), the trace of the stress-energy tensor (42), the effective action (which is \I/) 
(40), and the current (41b). With the help of the identity (B3) we see that this is indeed the same as T times the 



entropy from Eq. ( 73 1 . We have thus identified all terms in the generalized thermodynamic relation microscopically. 



B. Entrainment and superfluid density from field theory 



The independent degrees of freedom of our microscopic calculation are the chemical potential /i = d°tp, the tem- 
perature T = 0°, and the superfluid three- velocity v s — —Vip/fi, all measured in the normal-fluid rest frame, where 
the entropy current vanishes by definition, s l — 0. We have thus given 8 independent components out of the 16 
components of the 4 four- vectors j 1 *, s M , M , d^tp of the two-fluid formalism. The other 8 components are j M , s°, and 
0\ For the Noether current j M and the entropy s° we have field-theoretic and thermodynamic definitions. It remains 
to compute the spatial components of the thermal four- vector M . Additionally, we have to compute the coefficients 
A, B, C. They are defined as the derivatives of ^> with respect to the Lorentz scalars er 2 , O 2 , and O • dtp. However, 
this is not the form in which our is gi ven . Therefore, we need to find a different way t o compute these coefficients. 
This can be done with the help of Eqs. (58). First, we solve the spatial part of Eq. (58b) for Q l and insert the result 



into the spatial part of Eq. (|58a|). Together with the tem poral components, this yields three equations for the three 

where, for completeness, we also give the coefficients A, B, 



III 



variables A, B, C, whose solutions are listed in Table 

C, which are obtained from the inverse of Eqs. (60). With these results we immediately find 



A Vtp 



f + 



(vv-) 2 



(75) 
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Since all results are expressed in terms of quantities accessible from our microscopic calculation, we can for instance 
compute (for m = and in the low-temperature limit) 



A 



B 



An 2 T 3 1 - v 2 



16tt 4 T 5 1 - v 2 



W3/x (1 - 3v 2 ) 2 315^3^ (I - 3v2)5 



(25 + 78v 2 - 27v 4 ) . 



/r „ . 2 , 4tt 2 T 



l-v 2 8tt 4 T 6 65 + 256V 2 . - 402v 4 + 81 vf 

X { ~ VJ ~ 15y/3f (1 - 3v 2 )3 + 315^3^ ~ (1 - 3v 2 )6 



2^ 2 T 2 l-v? 



15^3 1 - 3v 2 315\/3> 2 



8tt 4 T 4 5 - 59v 2 + 27v 4 



27v° 



and 



~ - 



2v s r 

l-3v? 



1 - 



(l-3v 2 ) 4 



(l-v 2 )(5 + 9v 2 ) fnT 



7(1 - 3v 2 ) 



(76a) 
(76b) 
(76c) 

(77) 



We see that the coefficients A and C vanish at T — 0. This is in accordance to our zero-temperature discussion, where 
only B was nonzero. 



The connection between the coefficients A, B and the number densities n s , n n is given in Eq. (62). We can insert 



A and B as functions of n s , n n into the temporal component of Eq. (56b) to get also C as a function of n s , n n . The 
result is the useful translation 
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As a check, we see that B given in Table III is indeed the same as B — a/n s with n s from Eq. (50 1. It is now of course 



straightforward to also express A, B, C in terms of n s and n n . 

In the small-temperature approximation, the superfluid and normal number densities, measured in the normal-fluid 
rest frame, become 
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(79b) 



(Remember that n s is the superfluid density is the superfluid rest frame; the Lorentz factor n/cr = 1/vT — v s 
transforms it to the normal-fluid rest frame.) One can check that the sum of both densities gives the total charge 
density j° from Table |lj As expected, the normal density vanishes for T = and begins to increase with increasing 
temperature, while the superfluid density decreases. In a more complete treatment, we would expect the superfluid 
density to vanish at the critical temperature because the condensate melts. Remember that this melting is, in our 
one-loop effective action, a higher-order effect in the coupling constant, which we have neglected. The decrease of n s 
is therefore only due to the interaction between the two fluids. As an aside, note that for v s = the T 4 contributions 



in superfluid and normal densities cancel each other exactly. We have made this observation already below Eq. ( 46 ) 
where we have seen that in the m = limit there is no T 4 contribution to the density. 

Finally, we may express the generalized pressure W in terms of Lorentz scalars. This reformulation is instructive 
because it gives 'J in the form that is usually assumed in the two-fluid formalism. However, for field-theoretic 
calculations that start from an underlying Lagrangian or an effective action, this form of ^ can only be found a 
posteriori and is thus not of particular practical use. Our quantities in the normal-fluid rest frame T, /i, w s are 
translated into the relevant Lorentz scalars a 2 , O 2 , dip • via 
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We solve these equations for T, /i, and v s and insert the result into the effective action (47). Then, up to fourth order 
in the temperature we can write the generalized pressure as 



*(CT 2 ,e 2 ,e-<9V) 



4A 



90^ 



e- 



(81) 



■2v»v v , see Eq. (8.9) 



The term in the square brackets can be written as Q^Q^Q^ with the "sonic metric" Q 

of Ref . |11) . In other words, the Lorentz invariant T 4 term of the pressure in the presence of a superflow is obtained 
by replacing T 2 —> Q^Q^Q U in the T 4 term of the pressure in the absence of a superflow. In principle, we can use 



the higher order terms in Eqs. (80b), (80c) to write the T 6 contribution in terms of Lorentz scalars. However, we 
have not found a compact way of writing this contribution. But, we have checked that it is not simply given by the 
same replacement as for the T 4 term. This is no surprise since the sonic metric is constructed solely for systems for 
a Goldstone mode with linear dispersion. The T 6 term however, knows about the cubic term in the dispersion. 



VI. SOUND VELOCITIES 



As an application of our results we compute the two sound velocities of the superfluid. Related calculations can be 
found in the recent literature in the nonrelativistic context of superfluid atomic gases |48H51| . where the experimental 
observation of both sound modes is in principle possible, although challenging [52, 53 . Our results will give the sound 
velocities in the presence of an arbitrary superflow, i.e., an arbitrary relative velocity between the superfluid and the 
normal fluid (limited by a critical velocity, as our results will show). In particular, they will depend on the angle 
between the direction of the sound wave and the direction of the superflow. A similar calculation in the nonrelativistic 
context of superfluid helium has been performed in Ref. |54j where, in contrast to our calculation, the sound velocities 
are computed in the superfluid rest frame and without temperature corrections. We do find temperature corrections 
to the velocity of second sound, which, as we shall see, arise from the cubic terms in momentum of the dispersion 
relation of the Goldstone mode. 

The sound wave equations are derived from the hydrodynamic equations. We start from Eqs. (64), i.e., from the 
current and entropy conservation and the vorticity equation. [Of course, equivalently, one can start from the current 
conservation plus energy-momentum conservation ([9]).] In addition, we need the expression for from Eq. (55), 



which will allow us to rewrite derivatives of thermodynamic quantities in terms of derivatives of the independent 
variables. These are the chemical potential a = d°ip, the superfluid three- velocity v s (more precisely, V-0 = — A*v s ), 
the temperature T — 0°, and the normal fluid three- velocity v„. All these variables are now allowed to exhibit small 
oscillations in space and time about their equilibrium values, T — > T + <5T(x,t), a — > u + 5/x(x, t). We perform the 
calculation in the normal-fluid rest frame from the previous sections, i.e., the superfluid velocity has a (large) static 
and homogeneous equilibrium value on top of which the sound wave oscillations occur, v s — > v s + <5v s (x, i), while the 
static and homogeneous part of the normal velocity can be set to zero, v„ — > Sv n (x,t). Of course, we need to keep 
the oscillations of the normal fluid velocity <5v n (x, t) because there is no global rest frame in which they vanish. 

We employ the linear approximation in the oscillations. In this case, the temporal component of the vorticity 
equation is trivially fulfilled. From the remaining equations one can eliminate the normal velocity, such that one is 
left with two equations where the sound wave oscillations are solely expressed in terms of oscillations in T and fj, 
(oscillations of the superfluid velocity \/ip can be expressed in terms of oscillations of /x by applying a time derivative 
to the whole equation and using d$V^> — V/x). The derivation of the wave equations is quite lengthy, and we explain 
the details in appendix |D 1| One obtains the following system of two equations, 
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where u> = /in„ + sT is the enthalpy density of the normal fluid. Each term is a product of a space-time derivative - 
in which we can replace T by 5T(x, t) and /z by 5/x(x, t) - and a prefactor that only contains the equilibrium values 
T, fj,, and Vtp = -/iv s . 

Before coming to the general result, let us write down the wave equations in two limit cases. Firstly, let us set 
T = 0. In this case, the normal number density vanishes, n n — 0, and thus n — n s ji/a. With this relation and the 
zero-temperature expression n s = ct 3 /A (we set m = for simplicity in this subsection) one finds that all terms on 
the right-hand side of Eq. ( 82a ) vanish, and Eq. ( 82b ) can be compactly written as 



(83) 



With 



This wave 



Again, we recover the sonic metric — + 2v^v v , see remark below the generalized pressure (81 1 
5(i = 5fi e lk - x we obtain Q^k^ku = 0, which is Eq. (4.12) of Ref. [TT], see also Eq. (29) of Ref. [TB] 
equation has one physical solution uj = ui|k|, with the velocity of first sound u\. The solution is given in Eq. (85a 
(as we shall see below, this solution is unaltered by temperature effects up to the order we are working). 

Se cond ly, we discuss the limit case without superflow, Vtp = 0. In this case, only the first lines of Eqs. ( |82a 

with 8{i = (5yU e 1 (" Jt ~ k ' x ) and ST = <5T e l ^ t_k x ) we obtain two equations for the 



and (82bj) are nonvanishing. Now, 

Since we are interested in nontrivial solutions, we need to require the determinant of the 

n s + n n (which is true for Vip = 0) and 



two amplitudes Sfio, 5Tq. 

coefficient matrix to vanish. After a bit of algebra, using n 
resulting equation can be written as 
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This result is in exact agreement with the one given in Eqs. (19) - (22) of Ref. [43]. Now there are two physical 
solutions, uj = ii! 2 |k|, with the two sound velocities u\, 1*2 ■ The reason for the appearance of the second mode is that 
the presence of the second fluid component allows for relative oscillations between the two fluids. The solution of Eq. 



(84) is the v s = limit of the full result (85). 



In general, the full wave equations (82) yield very complicated results for the sound velocities. However, in our 
approximation for small temperatures up to order T 6 in the pressure, one can show that the resulting quartic equation 
for lu factorizes into two quadratic equations. This is explained in detail in appendix D 2 In this appendix we also 



explain that our truncation of the low-temperature series does not allow us to compute temperature corrections to 
the sound velo cities of order T and higher. The T 2 corrections, however, can be reliably determined. As one can 



see from Eqs. (D19), this is possible because of the T 6 terms in the pressure which originate from the |k| 3 term in 



the dispersion of the Goldstone mode. It turns out that there is a T 2 correction only to the second sound ui- 
explicit results are 
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where 9 is the angle between v s and the direction of the sound wave given by the wave vector k. As a consistency 
check, we confirm that u± is the (m = limit of the) linear part of the dispersion of the Goldstone mode from 
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FIG. 1: Velocities of first and second sound u\, U2 from Eqs. (85 I for six different values of the supcrfluid velocity |v s | between 
and l/\/3- All velocities are measured in the normal-fluid rest frame. In these polar plots, the sound velocities for a given 
angle between the direction of the wave vector and the superflow are given by the radial distance of the curve to the origin; 
the direction of the superflow is parallel to the horizontal axis and points to the right; the scale is normalized to the velocity of 
first sound in the absence of a superflow, as one can see in the upper left panel. The speed of first sound does not depend on 
temperature within our approximation. The speed of second sound is shown for three different temperatures: T — [(black) 
solid], T/n = 0.02 [(blue) dashed], T/ fi = 0.04 [(red) dashed-dotted] . For large superfluid velocities the temperature expansion 
breaks down, and we have only shown the results for T — 0. 



Eq. (32a I, which was computed as one of the poles of the propagator. In our approximation, this dispersion does not 
depend on temperature. We know that in general the dispersion of the Goldstone mode does become temperature 
dependent, see for instance Ref. |35j . where the melting of the condensate has been taken into account, however 
without any superflow. 

The velocity of second sound 1*2 becomes complex for certain angles 9 as soon as |v s | > 1/V3. Moreover, the T 2 
term of is divergent as |v s | approaches We have seen in Table [i] that all components of the stress-energy 

tensor and the current exhibit this divergence too. The expressions in that table show that due to this divergence 
the T 6 term (say, in the energy density T 00 ) becomes comparable to or even larger than the T 4 term for superfluid 
velocities close to (and below) the critical velocity 1/V3, even if T is very small. This suggests that a calculation 
to all orders in T must be performed to predict reliably the behavior in this close-to-critical regime. In analogy, for 
the speed of second sound close to the critical velocity and at nonzero temperatures we would also need a resummed 
result, and we cannot trust the truncated expression (85b I. 

We plot the two sound velocities for all angles and for various superfluid velocities in Fig. [I] Because of the 
breakdown of the temperature expansion we have just explained, the results for nonzero temperature are only shown 
up to a superfluid velocity where the T 2 correction is still smaller than the T = term. We see that both sound 
velocities are increased when they propagate parallel to the superflow and decreased when they propagate in the 
opposite direction. At T = 0, where the result can be taken seriously for all |v s | < l/VS, the speed of second 
sound decreases significantly when the critical velocity is approached, and goes to zero for all "backward" angles 
ir/2 < 9 < 3ir/2 (while the velocity of first sound only goes to zero for propagation exactly antiparallel to the 
superflow, 9 — tt). Interestingly, for a given superfluid velocity, the temperature effect always increases the speed of 
second sound for all angles. We know that for larger temperatures it must decrease again, because it has to vanish at 
the critical temperature where there is only one fluid in the system. Within our low-temperature approximation we 
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cannot see this decrease. 



VII. SUMMARY AND OUTLOOK 



We have discussed the dissipationless hydrodynamics of a relativistic superfiuid, starting from a complex scalar field. 
Our main goal has been to relate the field theory with the covariant two-fluid framework of superfluidity, motivated 
by quark superfluidity in the CFL phase and its astrophysical relevance. Our results can be summarized as follows. 

• Microscopic calculation. We have started from Bose-Einstein condensation in a (p 4 theory. The condensate has 
been assumed to be static and homogeneous (which corresponds to the simple hydrodynamic scenario of a static, 
homogeneous superflow). A crucial role is played by the phase of the condensate. While small oscillations of the 
phase correspond to the excitations of the Goldstone mode, rotations of the phase around the full U(l) circle 
give rise to a chemical potential (speed of the rotation) and a superfiuid three-velocity (number of rotations 
per unit length). In the presence of a non-zero superflow, the excitations of the Goldstone mode become 
anisotropic, and we have computed the resulting components of the conserved current and stress-energy tensor 
for nonzero temperatures. We have restricted ourselves to a weak-coupling, small-temperature approximation, 
including terms up to sixth order in the temperature. This has allowed us to present the microscopic results 
in an analytical form. For a study of arbitrary temperatures up to the critical temperature, a more elaborate, 
self-consistent calculation must be performed numerically. 

• Two-fluid formalism. In the relativistic two-fluid formalism, the basic variables are the charge current and the 
entropy current. In the dissipationless case, both are conserved. For both currents, we can define conjugate 
momenta. Now, if the two fluids are interacting with each other, neither of the currents is (four-)parallel 
to its own momentum, but also receives a contribution from the other momentum. This contribution is called 
entrainment, and it must be computed from an underlying microscopic theory. We have explained this formalism 
in detail, and in particular have shown the connection to other, equivalent, formulations in the literature. For 
instance, the entrainment coefficient (and related coefficients) can be expressed in terms of the superfiuid and 
normal-fluid charge densities (whose corresponding superfiuid and normal-fluid four-currents are, even in the 
dissipationless case, not separately conserved). 

• Relationship between them. There are several concepts and quantities in the two-fluid formalism that are 
usually not used in field theory, such as the generalized pressure that depends on Lorentz scalars. Therefore, 
one important aspect of this work has been a translation of these quantities in field-theoretic terms. For instance, 
we have proven that, once we assume that the microscopic calculation is performed in the rest frame of the 
normal fluid, it follows that the generalized pressure is given by the effective action. We have also expressed 
the coefficients that relate the currents with the conjugate momenta in terms of quantities that are well defined 
in field theory, namely the space-time derivative of the phase of the condensate, the components of Noether 
current, and the entropy. As a result, we have been able to compute these coefficients explicitly as a function 
of temperature, chemical potential, and superfiuid velocity. As an application, we have computed the angular 
and temperature dependent velocities of first and second sound in the presence of a superflow. To obtain any 
temperature correction it has turned out to be crucial to go beyond the linear approximation of the dispersion 
of the Goldstone mode. Cubic corrections in the dispersion give rise to T 6 corrections in the pressure and T 2 
corrections in the velocity of second sound, while the velocity of first sound remains temperature independent. 
We have found that the velocity of second sound increases with small temperatures although eventually, beyond 
our small-temperature approximation, it must go to zero when the temperature approaches the critical point. 

Our study opens up several directions for future studies. The most direct extension of our work is to go to higher 



temperatures, for instance in the CJT formalism, as we have briefly laid out at the beginning of Sec. Ill A In a 
more phenomenological context, our work can be considered as a first step towards a more complete picture of the 
hydrodynamics of color-flavor locked quark matter. There are several steps in this direction for the future. It would 
be interesting to include a (small) term in the Lagrangian that breaks the global U(l) symmetry explicitly and see in 
which sense superfluidity and the two-fluid picture survive. In a sense, this is similar to including dissipative effects 
because dissipation leads to a non-conservation of the entropy current, while explicit symmetry breaking leads to 
non- conservation of the charge current. It is also not unlike anomalous hydrodynamics, where the axial anomaly 
leads to a non-conserved current and which has been discussed in the recent literature 55-57], also in the context 
of superfluidity [.58-60 . The case of a broken symmetry, while also of theoretical interest, is relevant for the kaon- 
condensed CFL phase, a viable candidate phase of dense quark matter inside a compact star. In this phase the kaon 
condensate does break a symmetry that is, due to the weak interactions, not exact. It would also be interesting 
to include both superfiuid components of the CFL-A' phase, the one coming from kaon condensation and the one 
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coming from the spontaneous breaking of the U(l) associated with baryon number symmetry. It might be interesting 
to start from a fermionic formalism here, since baryon number is spontaneously broken by the quark Cooper pair 
condensate. Finally, our work can be taken as a starting point to include dissipation in the two-fluid picture and 
possibly improve on existing calculations of CFL transport coefficients. 
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Appendix A: Matsubara sum with anisotropic excitation energies 



Here we derive the result (30) for the effective action. The calculation shown here is formulated in a general 
way, such that it is also applicable to the stress-energy tensor and the current in Sec. |IIIB| In order to perform the 
Matsubara sum in Eq. (29), we write the determinant of the inverse propagator in terms of its zeros, 
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In the presence of a superflow Vi/>, the zeros are very complicated. The reason is the linear term in ko in the off- 
diagonal elements. With the help of Mathematica we obtain analytical, but very lengthy expressions for ef k . For 
small momenta, we can write 
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We now use the Matsubara sum 
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with j = 2 if i = 1 and vice versa, and an arbitrary function F(ko,\i) (without poles in the complex fco plane). For 
the effective action, F(k ,'k) is given by Eq. (31), for the stress-energy tensor and the current see Eqs. (44). Since 



we are interested in small temperatures, we may neglect the contribution from the massive mode, i.e., the two of the 
four terms in the sum where i = 2. More precisely: later, after writing coth[e^ k /(2T)] = 1 + 2/(e| k ) with the Bose 
distribution /, we shall only keep the thermal contribution, which, in the case of the massive mode, is suppressed for 
small temperatures. For the non-thermal (divergent) contribution, all terms have to be kept in principle. However, 
after renormalization, the contribution is subleading since it contains an additional factor of the coupling constant A 
and we shall neglect it. Therefore, after taking the thermodynamic limit, we can write 
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where, in the last step, we have changed the integration variable of the second integral k — s- 
F(k , k) = F(—k , — k) as well as 



-k, and have used 



(A6) 



This relation is easily checked for the small- momentum expressions (A2), and also holds for the full results. Due to 
the symmetries with respect to reflection of k, the poles e~ k have thus dropped out of the result, and the only physical 

excitations are e+ k . Therefore, in the main text, we have simply denoted e^.k = e+ k and £(k) = C + (k). 



Appendix B: Renormalization and useful identities for stress-energy tensor 



With the function from Eq. (38) we can write the effective action (27) as 
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On the other hand, using Eq. (29), we have 
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with the four- vector w M = (1,0,0,0) and A kl B k , C k given in Eq. (39). A useful relation between A k , B k , C k can be 
derived with the help of the explicit form of the determinant of the inverse tree-level propagator, 
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Next, we rewrite the stress-energy tensor. With 
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we can write the stress-energy tensor from Eq. ( 36a I as 
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where we have used the definition of ^ k (38) and have added a constant, diagonal tensor Y^ v which has to be 



doip = M- I n this case, because of the first line of Eq. (|B2|), the pressure P = yT becomes 



determined such that the conditions (37) are fulfilled. In o rder to implement these conditions we now set Vtp — and 
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where we have used the form of the pressure (Bl) and dko/dT = ko/T (due to the linear temperature-dependence of 
the Matsubara frequencies). Consequently, 
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On the other hand, the nonzero components of the stress-energy tensor without superflow are, from Eq. (B5|, 
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By comparing Eq. (B9) with (B6) and (BIO) with (B8) we conclude that = diag(l, 1, 1, 1). Inserting this into 
Eq. (B5), we can write the renormalized stress-energy tensor as 

/ f)TT \ 7 1 

T» v = -[ 2- g^U + - \Ckk^k v + B h d^d u ij> + A k (k^d u iP + fc"0"V) + ^u v \ , (Bll) 



with = (1,0,0,0). 



Appendix C: Small-temperature expansion 



Here we explain the small-temperature expansion needed in Sec. Ill C for the effective action, the stress-energy 
tensor, and the current density. We focus on the effective action in this appendix, but the other results are obtained 
analogously. 

Expanding in powers of the temperature corresponds to expanding the integra nd i n powers of |k|. In order to obtain 
the result up to T 6 , we expand the integrand of the momentum integral in Eq. (30) as 



(ei,k + ei ,-k)(ei,k + £2-k)( e i,k - e 2,k) 
and, for the dispersion in the argument of the Bose distribution, 



oilkl 



£i,k - ci|k| 



C2 



(CI) 



(C2) 



where at, 0,2, c±, C2 are angular-dependent, dimensionless coefficients. Inserting these expansions, introducing a 
dimcnsionless integration variable y = ci|k|/T, expanding in T/a, and performing the resulting integration over y 
yields 



where we have used the integrals 



dy 



e» - 1 



4A 



7T 

15 



2n 2 T 4 
45 



dy 



47T 



40?r 2 /a 2 2a lC2 \ T 
~7cf \T c7 



e« - 1 



8tt 6 
~63~ 



dy 



16tt 6 
~2A 



(C3) 



(C4) 



For the case without superflow, V?/> = 0, the angular integral becomes trivial. In this case, with dgip = fi, the full 
dispersions are given by Eq. ( 45 1 , and we have 



Cl 



I** 



C2 



\J H 2 — m 2 (3fi 2 — to 2 ) 5 / 2 



(C5) 
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and 



We thus find for the pressure 



ci 3ci 

ai = — , a 2 = — . (C6) 



T (/r^-m 2 ) 2 ^ 2 T 4 4c 2 7r 4 T 6 

p = — r ~ — — h 

V 4A 90c? 63/i 2 4 

{p 2 - m 2 ) 2 (3^ 2 - m 2 ) 3 / 2 7r 2 T 4 ^ 6 (3/i 2 - m 2 ) 1 / 2 4tt 4 T 6 
4A + ( M 2 -m 2 ) 3 / 2 90 (> 2 - m 2 ) 7 /2 63^ 2 



(C7) 



The expressions for the case with superflow are quite lengthy in general, and we give the final results in the limit 
in = in the main text, see Table [I] 



Appendix D: Calculation of sound velocities 
1. Derivation of sound wave equations 

We start from the hydrodynamic equations 

= d^, (Dla) 
= fya", (Dlb) 
= s^e" -d v ®»). (Die) 

Before we evaluate them, we collect some useful relations. We denote P = P n + P s = \& and thus can write with 
Eq. (55]) 

dP = .fdid^ + s^de^ 

n 7? ft) \ / in \ 

= ndfi + sdT - -± W • rfVV' + -y a • dSJip - s • d (^VV'J - s • d (^sj , (D2) 

where j — n, s° — s, d ip = 0° = T, and we have eliminated j and by using 

f = n n u» + n s ^ , 6" = -—d»i> + -u" , (D3) 
a s s 

where w = e n + P n = /in n + sT is the enthalpy density of the normal fluid. In the linear approximation, s times 
a space-time derivative is negligible, because s = s v n and we neglect products of v„ with space-time derivatives. 
Therefore, we may approximate 

Tl 

dP ~ ndp + sdT -W • dVip . (D4) 

a 



This relation is needed to express derivatives of any thermodynamic quantity in terms of derivatives of T, /i, and Vi/>- 

dP _ d ■ 



For instance, we can write don = dojf- = -§-8qP etc. and obtain the following useful identities, 



don = p do , + ^doT - . VA1 , (D5a) 

o\x op, a\x 

flb. = ^flbM+^flbT-^V^.V M) (D5b) 



- n l d(diip) o{diip) o{diip) a 

v '(» - -s^-g^+^w****^- (D5d) 



2G 



where doip — M has been used. With these preparations we can discuss the hydrodynamic equations. The current 
conservation (Dial obviously becomes 

~ don + n„V • v n — V • I — Vip 



(D6) 

where we have used u ~ v„ and V • [n n v n ) ~ n n V ■ v n . Inserting Eqs. ( D5a[ ) and (D5d) into this equation, taking 
the time derivative of the result, and multiplying the whole equation by u yields 



dn 



ds 



ji— dofi + (i—DqT — n s — Au + un n \7 ■ <9qV 



du 



du 



-M afa , )A dodifi + u ^^ )A d (^T - u ata /A djipdidjU . 



d(di1>) 



(D7) 



Due to the linear approximation, all expressions have the form (equilibrium quantity) x (second space-time derivative) , 
since products of two first space-time derivatives are of higher order. 
The entropy conservation (Dlbl reads 



~ d s + sV • v„ . 

Inserting Eq. (D5b|, taking the time derivative of the result and multiplying the whole equation by T yields 

c T^dla + T^dlT - T 9 ^!^ ■ Vd u + sTV ■ d Q v n . 



(D8) 



(D9) 



It is convenient for the following to add Eqs. (D7) and (D9), 
dn „ dr. 
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du ' T ~dT 
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9 V 



ds ds 
^d~u +T 



dT 



d$T- 



^ d{n s /<j) + T d{n s /<j) 



du 



dT 



<9(V?A) 



Vt/> ■ Vd u + 2/i ® S Vip ■ Vd T - 2^ ^"_ S (^ (V^ ■ V) 2 /i - n s -An , 



d{Vipy 



9(VV) 5 



where we have rewritten the derivative with respect to diip in terms of the derivative with respect to (V?/>) 2 . 
Finally, we need the vorticity equation (Die). The temporal component becomes 

= s • (VT + <9 ©) • 



(D10) 



(Dll) 



We can neglect this equation completely, since in both terms the normal fluid velocity is multiplied with a space-time 
derivative. The spatial components are 



~ s(d & + VT) ~ sVT + wd v n + sd a (^W>) 
where © from Eq. (D3| has been used. The last term needs some rearrangements, 
S do(^VV>) = 



(D12) 



s 



-Vif}dos + d Q (n n Vtp) = - — Vip d s + d (nVip) - d (n s -Vip) 

s \ a J 

-Vibdos + Vipdon + nVu — —Vipdou — udo ( — Vtp ) • 

a \ g ) 



(D13) 



Inserting Eqs. (D5a), (D5b), and (D5c| into this relation, the result into Eq. (D12|, and taking the divergence of the 

n n dn dn 



resulting equation, we arrive at 

~ n n Au + sAT + wV • dov r , 



n n ds ds 
~s~df ~ du~ 



2u 



ds 



• Vd u 



2u- 



dn 



s dT du ^d{Vip) 2 
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dT 



du 



2// 



■ Vd (i 

d(n s /a) 
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{Vip-Vyu. (D14) 



The normal fluid velocity can now be eliminated by solving this relation for V • c?ov„ and inserting the result into the 



other two equations: inserting it into Eq. (D9| yields, after multiplying the whole equation with w/(Ts), Eq. (82a 



in the main text, while inserting it into Eq. (D10) yields Eq. (82b). These are the sound wave equations from which 
the sound velocities are computed as follows. 
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2. Solution of sound wave equations 



Replacing the chemical potential and the temperature in all space-time derivatives of the sound wave equations 



Q by 6/j, = <5^ e iM ~ k ' x) and ST = ST e i( - ut -^ , the sound wave equations become 

~ [aiLj 2 + (a 2 + a^u 2 ^ 2 cos 2 9) + a 3 /x|v s |w cos 6] Su + (&i<I> 2 + b 2 + b 3 fj,\v s \u) cos 9) 5T Q , 

~ [A l u 2 + {A 2 



A 4 // 2 v 2 cos 2 9) + A 3 u\v s \Q cos 9] Sfx + (Biti 2 + B 2 + B 3 fi\v s \ui cos 9) 5T . 



(D15a) 
(D15b) 



where cos 9 = k • v s , Q = w/|k|, and we have abbreviated 



ai = — 



w dn 
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(D16a) 
(D16b) 



(D16c) 



(D16d) 



In general, the determinant of the coefficient matrix of the system of two equations (D15) yields a complicated quartic 
equation for G). However, we can simplify the result as follows. First one can check, for instance by explicit calculation, 
that the temperature dependence of the various coefficients is 



(4W 

a, = a; 1 



af ] T G 



Ai= A 



(0) 



A (4) T 4 



A (6) T 6 i 



& (3) T 3 



b^T 5 . 



B4 



B (3) T 3 



Bf ) T 5 , (D17) 



where i = 1, 2, 3, 4, j = 1, 2, 3, and where the prefactors in front of the various powers of T depend on the superfluid 
velocity v s . Since we have computed the pressure up to order T e , all terms of order T 7 and higher must be neglected 
in these expressions. We see in particular that only the A^s contribute in the limit T = 0. Now, for the determinant 
we encounter two kinds of products, namely Aibj and atBj, 



. 6 (5) T 2i 



T 3 L4f + A?>T* + A?>T«][b) , 
cuBj = T 7 [a 2 (4) + a^T^Fjf + Fiff 2 ] = 0(T 7 ) . 



T 3 A (<%(3) 



. 6 (5) r 2i 



0{T 7 ). 



(D18a) 
(D18b) 



The first line shows that the 0(T 7 ) terms are unknown in our expansion because a b^T 7 term in bj would give rise 

to a T 7 term in Aibj, but we have not computed b^p . Therefore, we must neglect all products of the form given in 
the second line since they are all of order T 7 and higher. In other words, it is consistent with our approximation to 
set ai ~ Bi ~ and use the T — results for Ai. In this case the quartic equation for u) factorizes into two quadratic 
equations, 



A^w 2 + {A\"> + A^'uWi cos 



,(o) 



(°)„2„2 



fj,\v s \Cj cos 6* , 







[ 6 (3) + b Wj*]Q* + [&W + ftWy*] + [ b 



(,(5)^,21-2 



(3) 



l (5)^21 



l(3) 



b { 3 5) T 2 ]u\v s \£bcosl 



(D19a) 
(D19b) 



After dividing out the overall factor T 3 of the quartic equation, the highest remaining power of temperature is 2, i.e., 
our approximation allows us to reliably compute the sound velocities up to T 2 . From the first equation we see that 
one of the solutions has no T 2 correction. This is the velocity of first sound. The second equation yields the velocity of 



28 



second sound which does have a T 2 correction. We see that the coefficients are needed to compute this correction. 
These coefficients arise from the T 5 terms in the entropy, i.e., from the T 6 terms in the pressure. Had we truncated 
our expansion of the pressure at order T , the velocity of second sound would have turned out to be independent 
of temperature. There are two physical solutions of Eqs. (D19), cj = Mi_2 |k| . For the explicit calculation of the two 
sound velocities u\, U2 we need n s and n n from Eqs. (791 and the entropy s (obtained by taking the derivative with 
respect to temperature of ^ = Tj_ from Table III). The result is given in Eqs. (851 in the main part of the paper. 
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